Exact results for intrinsic electronic transport in graphene 
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We present exact results for the electronic transport properties of graphene sheets connected to 
two metallic electrodes. Our results, obtained by transfer-matrix methods, are valid for all sheet 
widths and lengths. In the limit of large width-to-length ratio relevant to recent experiments, we 
find a Dirac-point conductivity of 2e^/\/3'i and a sub-Poissonian Fano factor of 2 — 3\/3/7r — 0.346 
for armchair graphene; for the zigzag geometry these are respectively and 1. Our results reflect 
essential effects from both the topology of graphene and the electronic structure of the leads, giving 
a complete microscopic understanding of the unique intrinsic transport in graphene. 
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Graphene, a graphite monolayer of carbon atoms form- 
ing a honeycomb lattice, has a distinctive electronic 
structure whose low energy excitations are described by 
masslcss Dirac fermions. The successful extraction of 
micron-scale graphene sheets from a natural graphite 
crystal, and their deposition onto an oxidized Si wafer 
was a truly seminal event which ushered in a new 
era of realistic experimental and theoretical exploration. 
The subsequent explosion of graphene activity has fo- 
cused on fundamental questions concerning the transport 
properties of relativistic particles in graphene and on its 
potential applications as a high-mobility semiconductor. 

Theoretical predictions Q for two-dimensional Dirac- 
fermion systems give an intrinsic conductivity ctq of or- 
der /h. Minimal conductivities around this value were 
observed at the Dirac point Q in Ref. while later 
measurements suggested that ctq 4e^/7r/i when the 
width-to-length ratio of the sample is sufficiently large. 
This is the value obtained using massless Dirac fermions 
and graphite leads in a Landauer-Biittiker (LB) formula- 
tion [5|, |6|- It is associated with a maximum of 1/3 in the 
Fano factor, Fq Q, which reflects the partial transmission 
of quantized charge through the finite graphene system. 
Measurements of the current shot noise in both ballistic 
and diffusive 0] graphene systems have indeed found 
that Fq « 1/3 in short and wide samples. Many authors 
have addressed different aspects of the graphene trans- 
port problem, which we summarize below. While the 
finite conductivity and suppressed Fano factor are gener- 
ally expected in graphene systems, the underlying physics 
remains rather poorly understood, not least because the 
carrier density at the Dirac point is zero. 

In this paper, we consider graphene sheets of both arm- 
chair (AGS) and zigzag (ZGS) geometry, connected to 
two metallic leads as illustrated in Fig. [1] By estab- 
lishing a transfer-matrix formulation within this minimal 
model, we present exact results for the anomalous intrin- 
sic transport properties of graphene. We demonstrate 
that the AGS and ZGS are completely different, and ex- 
plain in detail the non-universal dependence of a and F 




FIG. 1: Schematic representation of an armchair graphene 
sheet connected to electrodes at interfaces I and II. Primitive 
vectors a^; and for the electrodes and ai and a2 for the 
sheet give length L = VSN | ai | /4 and width W = (Af-l)|ai|; 
ai = 02 = ay with ai = 2.46A for graphene. The ZGS case, 
obtained by a 7r/2 rotation, has length L = (A'^ — 1)1^1 1 and 
width W = \/3A/|ai|/4, stiU with N x M sites. 



on geometry, filling, and gate voltage. 

The low-energy properties of graphene can be de- 
scribed by a nearest-neighbor, one-orbital tight-binding 
model for 7r-electrons on a hexagonal lattice. 
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where clj is an electron creation operator at lattice site 
Tij = {xi,yj), (...) denotes nearest-neighbor sites, ^2 is 
the hopping integral, and ^ the chemical potential. The 
two electrodes are represented by scmi-infinitc rectangu- 
lar strips with hopping to, while the interface hopping is 
ti . We take the interface contact to be perfect and impose 
open boundary conditions on the two free edges of the 
sheet; it is the geometry of these edges which determines 
our nomenclature (AGS or ZGS). Because graphene has 
two sublattices, sheets of size N x M lattice sites are 
taken to have width M and length N = Am (AGS) or 
= 2m (ZGS) with m an integer. 
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We begin with the AGS case (Fig. 1) by construct- 
ing a transfer-matrix equation for the scattering of elec- 
trons between two electrodes. In the Schrodinger equa- 
tion Hip{E) = Eip{E), the wave function is represented 

as tp{E) = = 4j|0) fo^' ^ij]^ '^i^h the 

complex cocfhcients to be determined. E ~ Ep is 
the Fermi energy of the electrodes, which is set by their 
occupation ric- There are M right- and M left-traveling 
waves (channels) in each electrode, each channel charac- 
terized by a transverse wavenumber ky ~ m'+i/2 '^ith 
n = 1,...,M. The longitudinal wavenumber fc" is re- 
lated to fc^ by Ep = -2to(cosfcJ + cosk^). 

With a unit-amplitude, right-traveling wave incident 
on the sheet in the nth channel of the left electrode, 

afj = E (<5„'„e'^"'"' +r„.„e-''="'-') sin(fc^'y,), 

" „' , (2) 

= Y.tn'ne'''- ^' sin(A;^ y^), 

n' 

for the left and right electrodes, where r„'„ and tn'n 
are respectively reflection and transmission coefficients 
from channel n to n' . For each site in the sheet, 
(S ~ P-'^ijy where T,5 specify the nearest- 
neighbor sites of Tij and jl = {jj, — Ep)/t2- We ex- 
press the M coefficients a for a given i as the vector 
di = {an, ttiA/)"^ in order to connect di with its neigh- 
boring slices through the 2Af x2Af transfer matrix 7i, 

(3) 

A translation period involves four different slices (Fig. 1), 
so % cycles through the four 2x2 block matrices 

A = 

c = \ ~^! 1, D=\~^\" -rn, (4) 

where X is lower-bidiagonal with nonzero elements equal 
to 1, y is the inverse of X, and / is the identity. 

Recursive application of Eq. ^ for N slices (A^/4 
translation periods) relates the coefficients and 
at the left and right interfaces through 



into the set of binary linear equations 
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where h = tl/{tot2),i2 = 



(5) 

By considering the 



one-period transfer matrix ABCD, one finds that the 
transverse modes in the AGS are unmixed by scattering 
processes, remaining independent and retaining the free- 
particle dispersion e„ = —2 — 2cosA:^. This makes the 
LB formalism underlying our transport calculations par- 
ticularly appropriate. Thus Eq. ([5]) can be decomposed 
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with a„ = (/i^ -/i'*)e„i 
c„ = (A^ - l)e,T^ and a„c„ -I- = 1. 

Analytic solution for t„„ from Eq. ([6]) gives the trans- 



mission probability T, 
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71 cosh (iV0„/2) + 72 sinh (iV0„/2) + 73 ' 



(7) 



where 71 = {v'l + + + 72 = sign(Kj_)(j/ii^2 + 

i'3'^4)/4, 73 = i-^f + A - 4 + f^l)/8, = (-2^-^ + 
^+i_i_cosfc")/sinA;", V2 ~ t _ cos fc^/ sin fc" , = 
Vi = C± = (K+/K- ± k_/k+)/2, k± = (c„ - a„ ± 
26„)i/2, f± = ii ± t2, and cosh6'„ = (a„ -f c„)/2. The 
conductivity a and Fano factor F may now be computed 
exactly from Eqs. dH) and ([7|), leading to 
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These expressions are completely general within the LB 
framework, and arc applicable for all sheet sizes {N, M). 

For the purposes of this Letter, we focus on the physi- 
cal insight contained in Eq. ([5]) for the situation relevant 
to most graphenc experiments, namely wide electrodes 
xitterncd onto the sample with rather narrow separation 
In this limit of large W/L{^ M/N), the sum in 
Eq. (j8|) is replaced by an integral over ky . For convenience 
wc set to = ti = t2, which creates no special symmetries. 
The Dirac-point conductivity ctq and the corresponding 
Fano factor Fq may then be expressed analytically as 

2A/3arctan(|cosfc^/sinfc^|) 
h TT sin fcJJ |cos fcj/ sin fcJl ' 
I sin fc^/cos fc^l 
2arctan ( | cos fcg / sin k^ \ ) 

where ky is the Dirac-point wavenumber 2tt/3 of the 
AGS and k^ is determined by Ep. At half- filling of the 
electrodes, i.e. Uc = 1 and Ep = 0, we obtain ctq = 
2e^/V3h sa 1.1547eV/i and Fq = 2 - SVS/tt ss 0.3460. 

Our results for the AGS are similar but not identical 
to the values Ae'^/nh and 1/3 of Ref. While the elec- 
tronic structure of the electrodes leads to a small quan- 
titative difference between the two studies, we will show 
below that the symmetry-breaking effect of the electrode 
interfaces causes a strong qualitative difference. The fact 
that tTo 7^ at the Dirac point p, = 0, despite the vanish- 
ing density of states, is an intrinsic property of the AGS 
quite distinct from conventional mesoscopic systems. 

To analyze the physical origin of this behavior. Fig. [5] 
shows the full dependence of ctq and Fq on ric and on the 
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aspect ratio M/N of the sheet for jl = 0. In Figs. 2(a) 
and (b), for ric = I, ctq and Fq alternate as M increases 
between semiconducting and metallic behavior, the lat- 
ter obtained when mod(2M + 1,3) = and there ex- 
ists a resonant channel with T„ = 1 at /cj^ = ky [in- 
set Fig. 2(a)]. As the sheet width is increased, the two 
branches merge at M/N ~ 1.5 with cto and Fq indepen- 
dent of M/N; only when M ^ N are sufficiently many 
channels with T„ ~ 1 available that their contributions 
to the sum in Eq. ([8]) are constant. For such sheets [in- 
set Fig. 2(b)], channels with T > 0.9 contribute to Fq 
with a distribution P{T) cx 1/Vl — T while channels 
with T < 0.1 have P(T) oc 1/T @. Ahhough P(T) 
resembles the universal bimodal distribution for a dis- 
ordered mesoscopic system, which also has Fq ^ 1/3, 
the underlying physics is completely different: the sub- 
Poissonian behavior is caused by the interference of rel- 
ativistic quantum particles, which results in transport 
contributions T„ ^ exp(— ]fc^ — ky\N) away from the res- 
onant channel [inset Fig. 2(a)]. This type of behavior, 
namely ctq 7^ and _fo ~ 1/3, is obtained in the AGS 
only for 0.63 < < 1.83 [Figs. 2(c) and (d)]. 

Figure 3 shows the effects of a gate voltage on a and 
F for ric ~ 1. For a finite sheet, the number of Fermi 
momenta (number of energy bands intersected) increases 
with fi^ each peak in a and F corresponding to one more 
contributing resonant channel. When the sheet is suf- 
ficiently wide [Fig. 3(a)], channels are added at nearly 
equal intervals, resulting in almost periodic oscillations, 
whereas for M/N < 1.5 [Fig. 3(b)] the effects of added 
channels appear quasi-periodic. Superposed on the os- 
cillation is a linear and slightly asymmetric behavior of 
a about the Dirac point. The former is a consequence 
of the linear dispersion of graphene and the latter of the 
electron- hole asymmetry caused by the electrodes [loj . 
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FIG. 2: (Color online) Dirac-point conductivity ctq (a,c) and 
Fano factor Fq (b,d) as functions of M/N for tic = 1 (a,b) 
and of ric with M/N = 10 (c,d). Dots denote ao = 2e^ /\/2>h 
in (c) and Fa — 2 — Sx/S/i" in (d) at ric = 1, while crosses 
denote cro — Ae^ /nh in (c) and Fo = 1/3 in (d), obtained 
at Uc ~ 1.39. Insets: transmission probability T{ky) around 
fc^ = 27r/3 in (a) and its distribution P{T) (see text) in (b). 
Calculations performed with A'^ = 1000 in the system of Fig. 1. 



Thus our exact results illustrate the inherent dependence 
of experimental observations on both W/ L and L [1, 0] , 
and demonstrate further that such behavior can be in- 
trinsic, rather than appearing only as a consequence of 
sample disorder or interfacial defects. 

We turn now to the ZGS. The geometry of this case 
requires a transfer matrix % expressed in terms of two 
2x2 block matrices and a quartic form of Eq. ^ for 
tnn, which is solved numerically to obtain a and F from 
Eq. ([8]). Figure 4 shows the dependence of cto and Fo on 
M/N, again with tg = ti = t2 and jl = 0. The ZGS also 
possesses metallic and semiconducting branches, which 
alternate with respect to the sheet length N rather than 
to its width M. The asymptotic behavior is metallic, 
with cr^ cx 16/73^ for mod(2iV + 1,3) = 0, and semi- 
conducting with o-g 0.2801/A^2 otherwise [Fig. 4(a)]. 
The corresponding Fano factors [Fig. 4(b)] are Fq ^ (1 — 
0.2801/iV2)/(l-H32.987VVf-) and F^ - 1 - 0.08223/A^^ 
the two branches merging only when ^ > 32.98N'^. 
Graphene sheets in this limit of W/L would therefore 
have (To = and Fq ^ I at the Dirac point, implying 
a Poissonian shot- noise quite different from the AGS. A 
finite minimal conductivity, do ^ 2e^/V3h (reaching its 
maximal value when Ep = ±1), and a sub-Poissonian Fq 
are obtained for all ric 7^ 0, 1, 2 [Figs. 4(c) and (d)]. 

The origin of the contrasting intrinsic transport prop- 
erties of the AGS and ZGS for the Dirac point lies in 
the special nature of zigzag chains in graphene. The key 
point is how this affects scattering at the interfaces. Be- 
cause g ai-f-Tj+s = for any site Vij in the sheet, the 
wavenumber of extended states is 27r/3 when projected 
onto the zigzag chain direction, and zero in the orthog- 
onal direction. In the AGS, zigzag chains are parallel to 
the interfaces so that ky = 27r/3 for mod(2M 4-1,3) = 0. 
Thus the incident traveling wave is not deformed at the 
interface and there is no interfacial scattering. Conse- 
quently, T{ky) = 1 and T{ky) ^ 1 in a regime of width 
0{1/L) about ky [inset Fig. 2(a)], resulting in a finite 
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FIG. 3: Conductivity a (a,c) and Fano factor F (b,d) for the 
AGS with M = 10000 and Uc — I, shown as functions of jl 
for M/N = 10 (a,b) and M/N = 1 (c,d). 
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FIG. 4: (Color online) cto (a,c) and Fo (b,d) for a ZGS, shown 
as functions of M/N at ric = 1 (a,b) and of ric for M/N — 
1200 (c,d). Inset: T{ky) around ky = 0. Red and blue curves 
indicate respectively metallic and semiconducting situations, 
calculated with iV = 1000 and 999. 



(To after the summation in Eq. ([S]) if 0.63 < Uc ^ 1.83 
[Fig. 2(c)]. In the ZGS, zigzag chains connect the left 
and right electrodes, ky = 0, and the armchair inter- 
faces involve two sublattices, with two values of fc" corre- 
sponding to each ky. This induces interfacial scattering. 
As a consequence, for Uc = 1 the transmission ampli- 
tudes are suppressed very strongly for any ky ^ ky and 
Tn = ^fej.fcs + [inset Fig. 4(a)], where is a very 
small background of width 0{1/L^) arising from interfa- 
cial scattering. Neither term contributes to the integral 
in the limit of large W/L and L, whence (To = and 
Fo = 1. When ^ 1 [18[, imaginary fc" values appear 
for some channels {ky}, giving contributions to T{ky) 
over a greater width and leading to a finite (Tq [Fig. 4(c)]. 
Thus it is the topological difference in the geometry along 
and across a hexagonal lattice which results in two fun- 
damentally different types of interfacial scattering, and 
hence in the contrasting intrinsic transport properties of 
AGS and ZGS systems. This microscopic insight was not 
included in any previous studies. 

Many investigations of graphene transport may be 
found in recent literature. Augmenting the general re- 
sults cited above, experimental studies of the conductiv- 
ity minimum have addressed the coherence of Dirac-point 
transport 11 1, the role of contacts and sample edges fl^ . 
and how interface charging leads to asymmetric gate- 
voltage effects |l3| . Many theoretical studies have consid- 
ered transmission coefficients in a finite graphene system, 
all restricted (as here) to the case of non-interacting elec- 
trons: from its weak interactions and the vanishing den- 
sity of states at the Dirac point, the fundamental trans- 
port properties of graphene are expected to emerge at the 
band-structure level. These investigations all differ from 
ours in the approximations applied, or in system size and 
geometry, or in the method of calculation, and hence in 
the nature of their conclusions. In an effective contact 



at the Dirac point makes all leads equivalent. Numerical 
treatments, of the same model fl5| and in a more gen- 
eral framework have probed size, gate- voltage, and 
impurity effects. While these and other studies 0, [l3| 
note that the AGS and ZGS cases should differ, the fun- 
damentally different nature ((To = 0, Fo = 1) of the ZGS 
case and the microscopic origin of the different intrinsic 
transport properties have been missed. Further, because 
we have analyzed the intrinsic transport arising due to 
lead and interface geometry, we may conclude that dis- 
order effects are not required to obtain the anomalous 
behavior observed in experiment 0, 0, [Hi . 

To conclude, we have presented exact solutions of the 
transfer-matrix equations for graphene sheets with metal- 
lic electrodes. Our results are microscopic and com- 
pletely general, and can be used to show that the Dirac- 
point conductivity and the Fano factor tend respectively 
to (To = 2e^/^/2>h and Fo = 2 — 'i^/Z/^T for armchair 
graphene sheets in the short and wide limit relevant to 
experiments. The same quantities tend to and 1 re- 
spectively for zigzag graphene sheets. Our exact results 
suggest that the measured finite minimum conductivity 
and sub-Poissonian Fano factor are the consequence of 
armchair rather than zigzag graphene systems, and show 
how this fundamental difference depends on the avail- 
ability of resonant transmission channels, which is deter- 
mined in turn by the geometry of the hexagonal lattice. 
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model [14| for a sufficiently large system, mode selection 
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